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Plants at high population density compete for light, showing a series of physiological responses known as the shade 
avoidance syndrome. These responses are controlled by the synthesis of the hormone auxin, which is regulated by two 
signals, an environmental one and an internal one. Considering that the auxin signal induces plant growth after a time 
lag, this work shows that plant growth can be modelled in terms of an energy-like function extremization, provided 
that the Markov property is not applied. The simulated height distributions are bimodal and right skewed, as in real 
community of plants. In the case of isolated plants, theoretical growth dynamics and speed correctly fit Arabidopsis 
thaliana experimental data reported in literature. Moreover, the growth dynamics of this model is shown to be consistent 
with the biomass production function of an independent model. These results suggest that memory effects play a non- 
negligible role in plant growth processes. 



1. Introduction 

Ecologists have been concerned with comp etition in 
plants since the beginning of the last century. 
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Clements 



19051) published the first formal definition of competition 
in plants, defining it as "the relation between plants occu- 
pying the same area and depen dent upon the same sup ply 



;upr. 

of physical factors" , and later (|Clements et all 119291 ) he 



considered that competition begins "when the immediate 
supply of a single factor necessary falls below the com- 
bined demand of the individual plant" . 

Since the 50s a number of authors have systematically 
studied the effe c t of c ompetition in plant systems. For ex- 
ample, iDonaldl (|195ll ) investigated in a seminal work the 
infra-spe cific light competitio n in pastures plants. More 
recently, IPacala et al. (1996) report experimental data 
providing little evidence of competition for nitrogen and 
water, and suggesting that the most important factor de- 
termining competition in sufficiently manured and watered 
plants is light. Accordingly, the present work will assume 
no stress due to insufficiency of nutrients or water (for 
models on dryland vegetation, see Kletter et al. ( 20091 ) and 
ref.s therein) and will study light competition in plants 
grown at high population density, i.e. experiencing den- 
sity stress. In particular, it will be concerned with intra- 
specific (inter-plant) light competition. The series of phys- 
iological responses demonstrated by plants experiencing 
density stress is collectively known as the shade avoidance 
syndrome and comprises effects such as stem elongation at 



expense of leaf and storage organ expan sion, inhib i tion o f 



branching, and acceleration of flowering (|Tao et all 120081 ). 

There are two main 'fingerprints' characterizing the 
height distribution of a community of plants under den- 
sity stress: 

i. bimodality: bigger individuals grow more than 
smaller ones, thus become bi gger and bigger, whilst s maller 
individuals grow much less (Ford and DiggleLll98lh . 

ii. hierarchy of exploitation (right skew): there is a 
large number of plants slightly smaller than the mean, and 
a small number of plants much greater than th e mean , 
i.e. the distribution is right (positively) skewed (|Harper[ 



1967). 



The present work is based on the consideration that 
growth responses in plants at high population density 
are controlled by the synthesis of several phytohormones , 
the most importa nt being auxin (jFriml and Saue 

E 120081: 

iTeale et all. 20061). Auxin ind uces plant growth after some 



delay dTaiz and Zeige 



suits (Tao et al 



200. 



2002) , and recen t expe rimental re- 



Stepanova et al. 1. 120081) show that 



the auxin pathway is mainly modulated by two signals, 
an environme ntal one and an inte rnal one (or, using the 
vocabulary of iGuedon et al.l ([20071 ) , ontogenetic) . The en- 
vironmental signal consists in the decrease in the red to 
far-red ratio of incoming light. In fact, plants grown at 
high population density perceive both absorption of red 
light by canopy leave s and reflection of far-red light from 
neighbouring plants ( Tao et al. . 20081) . In other words, a 
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plant 'recognizes' to be surrounded by other shading neigh- 
bours and responds producing auxin which induces growth 
after some time lag. On the other hand, the internal sig- 
nal is mediated by the hormone ethylene (jStepanova et al 
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2008), which stimulates the production of auxin even in 
isolated plants. 

A countless number of models have been developed to 
study both inter- and intra-specific competition among 
plants. Revie wing the diffe r ent sp atial community dy- 
namics models, iBolker et al.l (|2003l ) group them in three 
different frameworks, namely: 

a. Interactin g particle systems (IPS), which are stochas- 
tic, markovian ( Neuhauser . 200lh . continuous time models 
with discrete individuals located in t h e cells of a regular 
lattice (see, for example. |Diggle ( 19761) ; Durrett and Levin 
( 19941 ): iHendrv et al.l (Il996l )). 

b. Stochastic point processes, which track discrete 
individuals that compete locally, but assume that they 
occupy a single poin t in co n tinuous space (see, for ex - 
ample. iGandhi et all (Il998l) : IBolker and Pacalal (|l999l ): 
Dieckmann and Law! ( 2000 )). 

c. Patch models, which track group of individuals, al- 
lowing (noncontiguous patch models) or not (metapop- 
ulation or patch-occupancy models) multiple individu- 
als and multiple specie s per patch ( see, for ex ample, 
Levins and Culverl (Il97ll) : iGilad et al.l (12004120071) ). 

Strengths and weaknesses of each of the above frame- 
works have been exhaustively analyzed by IBolker et al 
( 20031) . 



In this paper, I propose a model that studies the old 
problem of the light competition in terms of the old 
method of the energy-like function extrcmization, but does 
not belong to any of the previous frameworks. In fact, even 
sharing some features with other IPS models, the model 
presented here is non-markovian. 

Indeed, markovianity is a very nice property in stochas- 
tic systems. It assures that in a system the future state 
depends onl y on the present st ate and is independent of 
past states ( van K amp en . 20071 ). In other words, the sys- 
tem has no memory and this considerably simplifies cal- 
culations. Unfortunately, auxin production induces plant 
growth after a time lag, and at the proper time scale this 
memory effect forbids modelling the system by exploiting 
the Markov property. 

I will show that in the present model the non- 
markovianity is necessary: 

• to reproduce plant height distributions that are bi- 
modal; 

• to correctly mimic experimental growth dynamics in 
isolated plants; 

• to be consistent with the biomass production model 
bv lCournede et all (|2008l) . 

I will consider the total amount of auxin present in a 
plant, but not its spatial distribution in plant's body, i.e. 
the well character iz ed ph enomenon of the auxin gradient 
(see Tanaka et al. ( 2006h and ref.s therein). For recent 
models on the auxin gradient the reader is referred to 
Jonsson et all (l2006l) ; [de Reuille et ail (|2006l ); ISmith et al ' 
(|2006l ): lNewell et al.1 (120081) . 



2. The model 

As mentioned in the introduction, auxin is the main reg- 
ulator of plant growth. For plants growing at high popu- 
lation density, auxin synthesis is regulated by: 

1) an environmental signal due to the absorption of red 
light by canopy lea ves and ref l ection of far-red light from 



2008) 



neighboring plants (|Tao et al. 

2) an internal signal, which makes plants grow- 
ing also when they d o not experience density stress 



(|Stepanova et all 120081) 



Consider a community of plants placed on a 2- 
dimensional N x N square lattice. The lattice has N 2 
nodes, each corresponding to a plant labelled as 
with i, j — 1, . . . , N. Each (i, j)-plant is associated to its 
own amount of auxin Sy(t) at time t. Each plant shades, 
i.e. interact with, its neighbours by which conversely it is 
shaded. I will consider only the interaction between first- 
nearest neighbouring plants. Indeed, trials with second- 
nearest neighbour interactions showed that the results do 
not change substantially. 

The auxin amount of a plan t is co nverted into biomass 
after a delay ( Taiz and Zeieer . 20021) that I will call a: if 
a certain auxin amount is produced by a plant at time t, 
then at time t + a the plant shoots to grow proportionally 
to the produced quantity of auxin (in this paper I will 
take into account only vertical shooting, i.e. I will assume 
the height of each plant to be proportional to its weight.) 
Accordingly, I will define the shade avoidance potential 
a,ij{t) as the function representing the predisposition at 
time t of the (i,j)-plant to avoid shade by shooting at 
time t + a. Its definition needs some steps. 

For simplicity, first consider an ideal subsystem of two 
plants, say the (i,j)- and the (i + 1, j)-plant, subjected 
only to the environmental signal 1 ). 

The predisposition at time t of the (i, j)-plant to shoot 
at time t + a depends first of all on its own auxin amount 
Sij(t). In fact, the greater sy(t) is, the more the 
plant has the possibility to grow at time t + a. 

Moreover, the (z,j)-plant perceives an environmental 
signal from the (i + 1, j)-plant, that absorbs red light and 
reflects far-red light. This signal "warns" the (i, j)-plant 
of the presence of its (i + 1, j) neighbour. The greater 
the (i + 1, j) auxin amount Si + ij(t) is at time t, the more 
(i + will shade the (i, j)-plant at time t + a. Because 
the (i,j)-plant "wants" to avoid the (i + shade, its 
predisposition at time t to shoot at time t + a has to de- 
pend also on the auxin amount Si+ij(t) of its neighbour. 
In fact, the greater Si+ij(t) is, the more the («, j)-plant has 
to grow to avoid the shade of the (i + 1, j)-plant at time 
t + a. 

Therefore, for a system of two plants subjected only 
to the environmental signal, the predisposition at time t 
of the (i, j)-plant to shoot at time t + a, i.e. the shade 
avoidance potential a,ij(t), reads 



Oij(t) = J Sij(t) ■ S i+ ij(t). 



(1) 



2 



J is a non-negative constant that expresses how strong 
plants shade each other or, in other words, how strong they 
have to compete for light. Because in a generic system 
of N 2 elements the (i,j)-plant has more than one first- 
nearest neighbour (namely, two if is at the corner 
of the lattice, three if it is at the border, four if inside), 
the (i, j)-plant shades and is shaded by all its first nearest 
neighbours. Therefore, all the contributions from the first- 
nearest neighbour have to be taken into account in the 
shade avoidance potential, i.e. 

aij (t) = J Sij {t) (2) 

n.n. 

where Sij(t) = s J _i J (i) + s l+ i 3 (t) + s. iJ _i(t) + s.y+i(i) 

is the sum over the first nearest neighbour of the 
plant. 

Eq. @ takes into account only the environmental signal 
1) perceived by the (i, j)-plant; if it had no neighbours, its 
shade avoidance potential would be zero, i.e. it would have 
no chance to produce biomass. Since isolated plants grow, 
the internal signal 2) can be taken into account by adding 
a term proportional to the auxin amount in the plant, 
hsij (t) , where h is a non- negative constant expressing the 
strength of the signal, i.e. how high the plant would grow if 
isolated. Therefore, the complete function expressing the 
shade avoidance potential of the generic (i,j)-plant reads 

a i:i (t) = h Sij (t) + J Sij (t) Sij (t). (3) 

n.n. 

Notice that the greater the shade avoidance potential is at 
time t, the more the (i,J)-plant will be able to shoot at 
time t + a, i.e. the more it will be able to intercept light 
and avoid shade. 

Studying the system in the mean field approximation, it 
is worth considering the global shade avoidance potential 
A(t) = J2iLiJ2f=i a ij (t) ■ Because all the plants of the 
community contribute with their mutual shading, A(t) is 
a global quantity that represents the predisposition of the 
whole system to avoid shade by shooting after the time lag 
or, in other words, the mean shade avoidance potential up 
to a N~ 2 factor. A(t) reads 

>4(*) = E(*«(*)[fc +J E *«(*)])> ( 4 ) 

ij n.n. ' 

where the sum J^ij = J2iLi Sj=i runs over au tne N 2 
plants. In passing, notice that Eq. Q is for mally analo- 
gous to the Ising hamiltoniar{3 (jHuangl . Il987t ). 

Accordingly to what discussed above, it is reasonable to 
assume that on average plants try to maximize the shade 
avoidance potential: the greater is the global potential 
at time t, the maximal will be the light exposure of the 
community at time t + a. 



1 Actually, it is the opposite. In fact, here the global shade avoid- 
ance potential will be maximized, while in the Ising model the hamil- 
tonian is minimized. 



Therefore, at each time step t the auxin amount of a 
randomly selected (i, j)-plant may be increased only if the 
global shade avoidance potential is maximized, as follows: 

Sy(*+1) = ««(*) (5) 
+ p£(t) e- bt Q[A(t)-A(t-a)]. 

£(t) G [— 1,1] is a randomly drawn number with flat 
distribution, p > the maximum updating amplitude 
for the selected auxin amount, 1/6 > the equilibrium 
timed. ©[■] is the Heaviside step function, defined as 
Q[x] = if x < and Q[x] = 1 if x > 0. The exponen- 
tial term e~ bt has been chosen to let the growing process 
dramatically slow down once reached the equilibrium time 
1/6. The positive discrete variable a represents the time 
lag introduced above and is calculated according to the 
flowchart in Fig. [T] Because at the very beginning of a 
plant life (dormant seed) the auxin amount is negligible, 
it is reasonable to choose as initial condition Sjj(0) = 0. 

Eq. ([5]) states that, at each time step, the selected auxin 
amount s^j can be updated by a randomly drawn quantity 
(which gets smaller and smaller as time goes by because of 
the exponential term) only if such an update implies the 
maximization (Q[A(t) — A(t — a)]) of the shade avoidance 
potential whose dynamics is defined as: 

A(t + 1) = A(t) (6) 
+ pm e- bt (h + jJ2sij(t)\ 

^ n.n. * 

with A(Q) = 0. 

£(t) is the same random number selected in ([5]). The sum 
^ runs over the nearest neighbours of the selected 
(i, j)-plant as above. 

A mortality rate has been introduced in the following 
way: if at certain time the (i, j)-plant has reached an auxin 
amount smaller than a threshold d, Sjj is set to zero and no 
longer updated during the simulation. This corresponds to 
the death of the (i, j) -plant. For the sake of clarity, the 
mortality process has not been shown in the flowchart in 
Fig.ffl ' 

Now, consider the function qij(t), that shall represent 
the height of the (i,j)-plant at time t. As discussed be- 
fore, the auxin amount of a plant is converted into biomass 
after a time lag a, and the biomass of a plant is assumed to 
be proportional to its height, qij{t). Therefore, setting the 
proportionality constant equal to one, before the equilib- 
rium is reached the plant height at time t corresponds to 
the auxin amount at time t — a. Then, at the equilibrium, 
the system does not evolve appreciably and, therefore, the 
plant height at time t corresponds to the auxin amount at 
the same time. In formulae, 

Hj(t) = s i:j (t-a) for t< 1/6, (7) 
?«(*) = Sijtt) for i> 1/6. (8) 



2 In this work, for equilibrium I mean the condition when the 
increments Asij become negligible. This happens at time t ~ 1/6. 
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Figure 1: Flowchart of the system dynamics. 



3. Numerical results 

I performed numerical simulations of the dynamics ([5][6| 
using the following parameter values: N = 10, p = .2, 
b = 1(T 4 , J = 10 2 , h = 1. The mortality threshold has 
been set to d = .3, which corresponds t o a rat e of m ortality 
of 70%, accordingly to experiments bv iFordl ( 19751 ). 

Fig. [U shows two height distributions at the equilibrium, 
i.e. at t > 1/6, for initial communities of iV 2 = 100 plants. 
Each histogram represents the height distribution of the 
dN 2 = 30 survived plants. 

In Fig. [2^, the simulation has been performed according 
to the model presented in the previous section: the sys- 
tem is non-markovian. The distribution res ults to be b i- 
modal and right-skewed (see experiments bv lFord ( 19751 )). 
Moreover, plants at the lattice border are on average (see 
Appendix |B|) shorter than central ones. 

On the contrary, in Fig. [2b, the time lag a has been 
set to be equal to zero, i.e. the system has been forced 
to be markovian. The height distribution looses bi- 
modality, showing that in this model framework the non- 
markovianity is a fundamental feature to interpret data of 
plants competing for light. 

The numerical simulations show that a varies with time. 
In particular, a can be consistently fitted by the simple 
saturating function 



a(t) 



m t 



(9) 



K + V 

where m is the plateau value reached by a and K = 1/ (26) 



a) non-Markovianity 




.2 .3 .4 .5 .6 .7 
% Heiqht 

b) Markovianity 
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% Heiqht 

Figure 2: Height distribution at the equilibrium, t > 1/b, for an 
initial community of N 2 = 100 plants. The rate of mortalit y have 
been set to 70%, accordingly to experiments bv IFordl J1975T) . The 
heights are reported in percentage relatively to the maximal height 
in each sample, a) Simulation of the system according to the model: 
the Markov property is relaxed. The distribution is bimodal and 
right-skewed, as in the experiments bv IFordl 1^975). b) Simulation 
of the system setting the time lag a = 0, i.e. forcing the system to 
be Markovian. The distribution looses bimodality. 



4. Isolated plants 

The model presented here is intended to mimic the be- 
haviour of plants at high population density. Indeed, it 
may be extended to isolated plants. 
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In Appendix \K\ I derived the expression for the growth 
(height) dynamics q(t) and the growth speed q(t) in iso- 
lated plants. They read 



«(*) 



P n b- 1 (1 -e 



p n 1 



TO if 



(if + t) 2 

where O is defined as in (|AB[) and 



-&/(*) 



/(f) = * - a(t) 



if -TO + i 

K/t+l 



(10) 
(11) 

(12) 



This allows to directly check the model consistency with 
data available in literature. 



Jouve et al.l ([19981 ) report experimental data for Ara- 



bidopsis thaliana floral stem elongation over time. They 
study circadian rhythms and report experimental data for 
Arabidopsis first inflorescence internode length and for 
first inflorescence extension rate as functions of time. The 
first inflorescence internode length can be assumed to be 
proportional to the plant growth dynamics q(t), and the 
first inflorescence extension rate, up to a constant, to the 
growth speed q(t). 

Fig.[3]shows q(t) (|10p fitted on Arabidopsis first inflores- 
cence internode length data and Fig. [4] shows q(t) ([TTjl fit- 
ted on Arabidopsis first inflorescence extension rate data. 

The agreement are satisfactory in both cases; fluctua- 
tions in first inflorescence extension rate data are due to 
the circadian rhythm, which is not taken into account by 
this model and which can be easily av eraged out. 

Interestingly, ICournede etaD ([2008) propose the follow- 
ing empirical formula for the biomass production Q(t), 



Q(t)=0\(t) (l- e -^W), 



(13) 



where /3 and 7 are characteristic parameters, S(t) the leaf 
surface area of the plants, and X(t) a function of environ- 
mental conditions related to the evapotranspiration that I 
here assume to be constant. 

Considering again the Arabidops i s case , experimental 



data from Fig. 3a in ICookson et al.1 (|2007j) show that, in 



the time window of interest (0-8 days), Arabidopsis leaf 
surface area over time can be consistently fitted (reduced 
X 2 ~ 1) by af(t), being a a proportionality constant. 

This suggests that the growth dynamics q(t) obtained 
here and the bioma ss production Q(t) proposed by 
Cournede etaD (12005' 



are consistent. While parameter 
identification between the two models needs further in- 
vestigation, it is worth observing that f(t) can correctly 
reproduce experimental data for leaf surface area over time 
only in the case m/0. In other terms, it is necessary the 
system to keep memory of the time lag a even in the case 
of isolated plants. 
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Figure 3: Arabidopsis thaliana growth dynamics. Eq. 111! (continu- 
ous line) and experimental data from Fig. 1, dataset A (12 hours of 
dark every 12 hours of light), in I Jouve et al.l d 1Q9ST> (dashed line). 
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Figure 4: Arabidopsis thaliana growth speed. E q. (] 1 [I fco n tinous 
line) and experimental data from Fig. 2A in I Jouve et al.l (l998) 
(dashed line). 



5. Sensitivity analysis 

Focusing back to the general model (J > 0, N > 1), it 
is interesting to analyze its parameter sensitivity. 
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I have performed a stochastic analysis of the model 
(Appendix IB]) . whose main result is, according to (|BH3|) 
and (TEH]) . 



(sy) S i ~ function(/o, 6, N). 



(14) 



This means that on average the auxin amount in a generic 
(«,j)-plant does not directly depend on J and h. Numer- 
ical simulations confirmed this result. Speculating, this 
suggests that plant height does not directly depend either 
on the R:FR ratio or on the hormone ethylene density. 
This observation is c onsistent with the fa ct that plant sen- 
sitivity b oth to light (jMaloof et al"l|200ll) and to hormone 
ethylene ( Chen and Bleecken , 19951) stimulations rapidly 
saturates. 

6. Discussion 

This paper proposes a simple model for the growth at 
high population density of plants competing for light. The 
model is based on the consideration that the plant growth 
is controlled by the auxin synthetic pathway, which is reg- 
ulated by two signals, an internal one related to the hor- 
mone ethylene, and an environmental one, due to the mu- 
tual shading. Moreover, the auxin signal is converted into 
biomass only after a time lag, a, that varies on time. This 
memory effect forbids modelling the system markovianly 
at the natural time scale. In fact, a single plant has the 
predisposition to avoid the shade experienced at time t 
only by shooting at time t + a. At the system level, this 
can be modelled by maximizing the global shade avoidance 
potential (gj). 

Consequently, this model cannot be considered to be- 
long to any of the three frameworks of spatial community 
models drafted bv lBolker et"aH (j2003h . namely the Inter- 
acting Particle Systems, the Stochastic Point Processes, 
and the Patch Models. In fact, even sharing some features 
with IPSs, the model proposed here is non-markovian. In- 
deed, it can be considered closer to those models that de- 
scribe in terms of energy-like function(al)s a broad range 
of life sc ience phenomena, such as X chromosome inac- 



tivation (Nicodemi and Prisco, 120071) . epithelial cell po- 
larization (IVeglio et all. 120091) . and eukaryotic chemotaxis 
(|Ferraro et all l2008t iGamba et al.Ll2009h . 

The present work shows that the simulated plant height 
distributions are bimodal and right skewed, as in real com- 
munities of plants at high population density. On the con- 
trary, if no delay is assumed between the auxin signal and 
the biomass production, i.e. if the time lag a is set to zero, 
simulations show that the distribution of plant heights is 
unimodal and cannot reproduce real distributions. 

Extending the model to isolated plants, an analytical 
expression is obtained both for the growth dynamics (1101) 
and the growth speed They have been fitted to Ara- 
bidopsis thaliana first inflorescence length data and on first 
inflorescence internode extension rate data, respectively. 
Theoretical functions and experimental data show a good 
agreement. 



Moreover, the obtained growth dynamics (1T0|) is con- 
sistent with the bioma ss production function proposed by 
Cournede et~al] (|2008l) . 



The bimodality in height distribution, the good match 
between experimental data and theoretical functions for 
isolated plants, and the consistency with an empirical 
model independently derived, can be considered the main 
results of this work. They hold as long as the time lag a is 
taken into account. If the model is forced to be markovian, 
a — 0, real distributions and dynamics can no longer be 
mimicked. This suggests that memory effects play a key 
role in the plant growth process and that they cannot be 
neglected. 

Interestingly. lHara ( 1984 ) studies the dynamics of plant 
populations with a stochastic model based on the Kol- 
mogorov forward equation, which needs markovianity, and 
observes that the Markov property applied to plant growth 
is reasonable as a first approximation and needs further 
investigation. In fact, it seems not to be a case that his 
model does not reproduce bimodality in plant height dis- 
tributions. 

A further step in the comprehension of plant growth 
would be to experimentally check the effects of varying 
the time lag between the auxin production and its conver- 



sion in biomass (for example, along the lines of iTao et al 



(2008)). both in high population density and isolated 
plants. Accordingly, from a theoretical point of view it 
would be interesting to study the bimodal-to-unimodal 
switch in height distribution due to the loss of markovian- 
ity, in parti cular in terms of statistical physics of ph ase 
transitions (jHuangl . 1987t iBarabasi and Stanley . 19951 ). 



Appendices 

A. Analysis of the Dynamics 

To analyze the dynamics, it is worth switching from dis- 
crete to continuous time by substituting t — > N 2 t. This 
yields ([5]) to read 



Sij(t) = p e- mH mj(t), 
where 

m 3 {t)^m B[A(t) - A(t - a)] 



is a random term. I will treat it as a gaussian non-zero- 
mean colored noise with time-independent mean, which I 
shall label as 

Ski ee (t)) = <£(t) 6 [A(t) - A(t -a)]). ©5) 

Therefore, the average of (|AHp over different system's re- 
alizations is simply 



(s t3 (t)) = P n ij e~ bN2 \ 

with the straightforward solution 



IB) 



which has the following steady state, 



= lim (sij(t)) = 



p 

b N 2 



423) 



The dynamics of an isolated plant height over time q(t) 
may be defined as 



q(t) = (s(t)) 



t=t-a(t) 



where (s(t)) is the same as in (jAttl) . but without position 
indexes because it is referred to an isolated plant. It is cal- 
culated at time t — a(t) to reflect that auxin is not imme- 
diately converted into biomass, as discussed in Section [2] 
(|AI7P is the same as ([7]) except for the average, which is 
intended to reduce fluctuations due to the isolated plant 
condition. The shade avoidance potential is ideally eval- 
uated in the condition J = 0, i.e. no interaction, and 
N= 1. 

Using ©, ([AT)]) and (LATj) . q(t) reads 



lit) = <«(*)> I «=t_„ (t) 



p n b- 1 (i 



and 



q(t) = p n (l 



m K 
(K + tf 



-b(t-mt/(K+t))\ 



-b(t-mt/(K+t)) 



The fact that a tends to the finite value m implies that 
at the equilibrium the delay gets negligible, letting q(t) = 
(s(t)) and q(t) = (s(t)) (|Aljl . This is consistent with 
(JHJ), that holds after the equilibrium. 

B. Analysis of the Parameters 

In order to explain the model dependence on the param- 
eters, it is worth evaluating analytically the noise mean re- 
introduced in (|AB|) . Its definition is general and valid for 
the whole system of plants. 

(jAll) is a Langevin-like equation with a non-zero-mean 
colored noise, no drift term and time dependent diffusion 
term. It is possible to perform some transformations to 
obtain a more standard Langevin-like equation. Defining 



i(t) = s i:j (t)e 



bN^t 



and inserting (IAH[) in the time derivative of (jBHI) , I obtain 
*«(*) =Fij(x)+p fj i:j (t), fBfc) 
where 



Fij(x) = bN 2 x l0 {t) +pQij, 

fjij(t) = Tjij(t) -Clij. 



IB) 



Consider that x(t) and s(t) have, at the steady state, the 
same pdf. 

(|BBI) is a Langevin-like equation with time independent 
drift term, constant diffusion term and zero-mean colored 
noise. 



Assu ming that such a no ise is a Ornstein-Uhlenbeck pro 
cess (|van Kampenl . 120071 ). i.e. 

dfjij (t) 



dt 



C(<), 



where /Lty is the correlation time and £(t) is a white noise 
with zero mean and delta time-correlation, it is possible to 
write the Fokker-Planck equation corresponding to (IBE 
i.e. 



dpij{x,7],t) 
dt 



= - bN 2 + 7p ij (x,r),t) + 



+ IV 



dpij(x,r),t) 1 d 2 pij(x,r),t) 
2 ' 



drj 2 drf 

Even at the steady state, (|BB|) is not easy to be analytically 
solved. However, for a huge set of parameters (|BB|) lies in 
the validity ra nge of the Unified Colored Noise Approxima- 
tion UCNA (Jung and Hanggi . 19871 lllanggi and Junel 
19951: lHasegawal l2008h. In f act to apply the UCNA it is 



necessary (jHanggi and ■Tund . ll99^ that 



> 



According to (|BH[) and to the fact that Sj 
(|BI7[) corresponds to 



1, constraint 



t > 



logjfiij/b) 
bN 2 



t tt: • 



Therefore, since equilibrium time 1/6 3> the UCNA 
can be applied for the parameter set used in the reported 
numerical simulations. Applying UCNA, the following ef- 
fective Langevin equation holds 

xn (t) ~ o-y (x) + p 4C(*) . (P>) 

where is again a white noise with zero mean and delta 
time-correlation, while I defined 

I calculated the correlation time /xy by averaging the nu- 
merical values of a over all the lattice nodes and over dif- 
ferent realizations of the system. As an example, labelling 
the site at the center of the lattice as cc, the sites at the 
center of the lattice borders as bb and the sites at the cor- 
ners as aa, I obtained, in continuous time units, fi cc ~ 54, 
p bb ~ 50 and fi aa ~ 47. 

In passing, notice that the drift term of (|BP[) . cr^Fy (cc), 
is linear in Xij\ therefore this model lies in a class of 
universality different from t he principal growth models^ 



(Barabas i and Stanley! . 119951 ). Then, it is possible to write 



the Fokker-Planck equation corresponding to (|B9[) , 



dpij(s,t) 
dt 



2 P ^ 



ds 

4 d 2 p l0 (s,t) 



8s 2 



©LI) 



^In particular, in the Kardar-Parisi- Zhang model iKardar et all 
1986) the drift term is equal to v\7 2 Xij + ^(Vxij) 2 , where v and /i 
are characteristic parameters of that model. 
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Figure 5: Eq. ((AQ (1 

incs) and numerical results from the simulation 
of the dynamics (symbols). Arbitrary units. 



and consequently to calculate the probability distribu- 
tion function at the steady state Pij(s) = lim^oo Pij(s,t) 
(jvan Kampenl . l2007t ) 



Pij (s) 



Therefore, 



( s ij) 



ij 1st 



e*p[-p- 2 a^(bN 2 s 2 + 2pO, 3 s)] 
/o 1 dy exp[-p-V:- 3 (WVV + 2p^y)\ 



ds s Pij(s) = g{Vtij,. . . ), 



4112) 



(Hl3) 



where the last term is to express that the analytical so- 
lution of the integral in (|BIL3[) . that I do not report for 
reasons of shortness, is a function of the noise mean Qij. 
Comparing (|BH3p to (|AB|) the following implicit equation 
in Vtij holds 



©4) 



which has solutions, for plants in different position of 
the lattice labelled as above, ft cc ~ .23, flbb — -20 and 
fl aa ~ .18. In Fig. [5] analytical solution (IA5|) and nu- 
merical results are plotted. 
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